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Abstract 

The changes of dynamical behaviour of a single fullerene molecule inside an armchair carbon nanotube caused by 
the structural Peierls transition in the nanotube are considered. The structures of the smallest C 20 and Fe@C 20 
fullerenes are computed using the spin-polarized density functional theory. Significant changes of the barriers for 
motion along the nanotube axis and rotation of these fullerenes inside the (8,8) nanotube are found at the Peierls 
transition. It is shown that the coefficients of translational and rotational diffusions of these fullerenes inside the 
nanotube change by several orders of magnitude. The possibility of inverse orientational melting, i.e. with a 
decrease of temperature, for the systems under consideration is predicted. 



Introduction 

The structure and elastic properties of carbon nanotubes 
are studied in connection with the perspectives of their 
applications in nanoelectronic and nanoelectromechani- 
cal devices and composite materials, and are also of fun- 
damental interest, particularly for physics of phase 
transitions. For example, superconductivity [1], com- 
mensurate-incommensurate phase transition in double- 
walled nanotubes [2], spontaneous symmetry breaking 
with formation of corrugations along nanotube axis [3] 
and structural Peierls transition in armchair nanotubes 
[4-9] have been considered. In the present Letter, we 
consider a fundamentally new phenomenon related to 
phase transitions in nanosystems. In other words, we 
consider the possibility of inverse orientational melting 
for molecules encapsulated inside nanotubes caused by 
structural Peierls transition in the nanotubes. 

The possibility of Peierls transition in carbon nano- 
tubes was first considered in [4]. As a result of this tran- 
sition, armchair nanotubes become semiconducting at 
low temperature, and Peierls distortions lead to the 
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Kekule structure (see Figure 1) with two essentially dif- 
ferent C-C bond lengths and a triple translational period 
(three times more hexagons in the translational unit 
cell). In previous studies, the Peierls gap [5-7] and the 
temperature of the transition to the metallic phase with 
equal C-C bond lengths [4,5,8] were estimated. Recently, 
the Kekule structure was calculated for the ground state 
of an infinite armchair (5,5) nanotube by PM3 semiem- 
pirical molecular orbital calculations [9]. It was shown 
that, for the (5,5) nanotube, the difference between C-C 
bond lengths for semiconducting phase is 0.03 A, 
whereas the difference between nonequivalent C-C bond 
lengths for metallic phase is only 0.006 A [9]. 

Note that density functional theory (DFT) calculations 
for the (5,5) nanotube of a finite length [10,11] also gave 
a 60 atom periodicity of physical properties on the 
length of nanotube segment which is consistent with 
Kekule structure for infinite armchair nanotubes. More- 
over, X-ray crystallographic analysis of chemically 
synthesized short (5,5) nanotubes [12] shows the Kekule 
bond length alternation pattern, which was in good 
agreement with DFT and PM3 calculations also per- 
formed in [12]. By the example of the infinite (5,5) 
nanotube, it was demonstrated that the structural Peierls 
transition connected with spontaneous symmetry break- 
ing takes place not only with an increase of temperature, 
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Calculated Kekule structure corresponding to the ground state of the (8,8) nanotube Nanotube axis is shown by the dashed 



Figure 1 

line. 



but also can be controlled by uniaxial deformation of 
armchair nanotubes [9]. 

A dynamical behaviour of molecules encapsulated 
inside nanotubes can correspond to the following 
regimes: oscillations about a fixed position and/or a 
fixed orientation of the molecule (regime A), hindered 
motion along the nanotube axis and/or rotation of the 
molecule (regime B) and free motion and/or rotation of 
the molecule (regime C). In the present Letter, we show 
the possibility of changes of the dynamical behaviour of 
molecules encapsulated inside armchair nanotubes as a 
result of the Peierls transition in the nanotube structure. 
In other words, these changes can include switching 
between the regimes A and B, switching between the 
regimes B and C, and the changes in diffusion coeffi- 
cients corresponding to the hindered motion and/or 
rotation of the molecule (regime B). The considered 
changes are possible in the case where the regime B 
takes place for at least one phase of the nanotube, i.e. 
the temperature T P of the Peierls transition should cor- 
respond to the temperature range of the regime B (the 
hindered motion and/or rotation of the molecule). In 
other words, the temperature T P should be of the same 
order of magnitude (or a few orders of magnitude less) 
as energy barriers AE for motion and/or rotation of the 
molecule inside the nanotube at this phase. Note that 
inverse melting of motion and/or rotation of the mole- 
cule is possible, if the Peierls transition from the high- 
to low-temperature phase of the nanotube occurs with 
switching from the regime B to the regime C or switch- 
ing from the regime A to the regime B. 

Estimations showed that the Peierls transition tem- 
perature is T P - 1-15 K [4,5,8]. According to calcula- 
tions [13,14], the barriers of the value close to this 
temperature range were obtained for rotation of the full- 
erene C 60 inside the C 60 @C 2 4o nanoparticle. It was also 



found that the changes of bond lengths of the fullerene 
C2401 the outer shell of these nanoparticles, within 0.06 
A lead to an increase of the barriers for rotation by 
more than an order of magnitude [13,14]. The changes 
of the nanotube bond lengths caused by the Peierls dis- 
tortions are of the same order of magnitude (about 0.03 
A for the (5,5) nanotube [9]). Note also that the size of 
an encapsulated molecule, and therefore, the nanotube 
radius cannot be too large, since the magnitude of the 
Peierls distortions decreases with an increase of the 
armchair nanotube radius [6]. 

Thus, taking into account the above considerations, 
we have chosen the smallest fullerene C 20 and the mag- 
netic endofullerene Fe@C 2 o to investigate changes in the 
dynamical behaviour of molecules inside nanotubes at 
the Peierls transition. It has been shown that the (8,8) 
nanotube is the smallest armchair carbon nanotube 
which can encapsulate the fullerene C 20 [15]. A carbon 
nanotube with the fullerene C 20 inside was also used as 
a model system to simulate a drug delivery via the 
nanotube [16]. 

This Letter is organized as follows: "Fullerene and 
nanotube structures" section presents the DFT calcula- 
tions of the structure of the C 20 and Fe@C 20 fullerenes 
and the PM3 calculations of the structure of the (8,8) 
nanotube. "Fullerene-nanotube interaction" section pre- 
sents the semiempirical calculations of the barriers for 
motion and rotation of the fullerenes inside the nano- 
tube. The section that succeeds the latter is devoted to 
the dynamical behaviour of molecules inside the nano- 
tubes. Our conclusions are summarized in the final 
section. 

Fullerene and nanotube structures 

Structures of the C 20 and Fe@C 20 fullerenes have been 
calculated using the spin-polarized density functional 
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theory implemented in NWChem 4.5 code [17] with the 
Becke-Lee-Yang-Parr exchange-correlation functional 
(B3LYP) [18,19]. Eighteen inner electrons of the iron 
atom are emulated with the help of the effective core 
potential - CRENBS ECP [20] (only 8 valence s-d elec- 
trons are taken into account explicitly). The 6-31G* 
basis set is used for describing electrons of the carbon 
atoms. 

The ground state of the fullerene C 20 is found to be a 
singlet state and has D 2 h symmetry. The calculated 
energy of the triplet state of the fullerene C 20 is found 
to be 64 meV greater than the energy of the ground 
state. The ground state of the endofullerene Fe@C2o is 
found to be a septet state and has C 2 h symmetry. 

The calculated structures of the ground states of the 
C 20 and Fe@C 20 fullerenes are shown in Figure 2. The 
iron atom locates in the centre of the endofullerene. 
The smallest and the greatest distances between the car- 
bon atoms and the C 2 o fullerene centre increase by 6 
and 1%, respectively, as a result of the iron atom 
encapsulation. 

The semiempirical method of molecular orbitals modi- 
fied for one-dimensional periodic structures [21] with 
PM3 parameterization [22] of the Hamiltonian has been 
used to calculate the structure of the (8,8) nanotube. 
The method was used previously for calculating the 
Kekule structure of the (5,5) nanotube ground state and 
for studying structural transitions controlled by uniaxial 
deformation of this nanotube [9]. The adequacy of the 
PM3 parameterization of the Hamiltonian has been 



demonstrated [23] by the calculation of bond lengths of 
the C 60 fullerene with I h symmetry: the calculated values 
of the bond lengths agree with the measured ones [24] 
at the level of experimental accuracy of 1(T 3 A. The cal- 
culated Kekule structure of the (8,8) nanotube ground 
state is shown in Figure 1. The difference between the 
lengths of short and long bonds of this Kekule structure 
of the (8,8) nanotube is close to such a difference of the 
(5,5) nanotube [9]. The Peierls distortions include also 
radial distortions of the armchair carbon nanotube with 
periodicity of half of the translational period of the 
nanotube (for details see [9]). In the case of the (8,8) 
nanotube, the longest nanotube radius is 0.547 nm, 
while the shortest radius is 0.544 nm. 

Fullerene-nanotube interaction 

The structures of the fullerenes C 20 and Fe@C 20 
obtained in "Fullerene and nanotube structures" section 
have been used for finding the ground state position 
and for studying motion and rotation of these fullerenes 
inside the (8,8) nanotube. Two structures of the (8,8) 
nanotube have been considered: the Kekule structure 
calculated in "Fullerene and nanotube structures" sec- 
tion, and the structure with all equal bonds 1.423 A in 
length (so that to be equal to the average bond length 
of the calculated the Kekule structure). The calculations 
of the interaction energy between walls of double-walled 
carbon nanotubes showed that the barriers for motion 
of the short wall relative to the long wall are very sensi- 
tive to the length of the long wall [25]. Thus, the size of 
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Figure 2 Calculated structure of the C 20 (a) and FeiaC 20 (b) fullerenes. The bond lengths are in angstroms. The distances between the 
carbon atoms and the fullerene centre (in angstroms) are denoted in italics. The atoms which have smaller and greater distances to the 
fullerene centre are shown by the open and filled circles, respectively. C 2 symmetry axis is shown by the dotted line. 



Poklonski ef al. Nanoscale Research Letters 201 1, 6:216 
http://www.nanoscalereslett.eom/content/6/1/216 



Page 4 of 10 



the system is too large for ab initio calculations. The 
analogous problem exists also for the considered case of 
a single fullerene inside a nanotube. Therefore, the 
interaction between carbon atoms of the fullerenes and 
the nanotube at the interatomic distance r is described 
by the Lennard-Jones 12-6 potential 

L/ w = 4 £ ((^) 12 -(^) 6 ) (1) 

with the parameters e = 2.755 meV, a = 3.452 A. These 
parameters of the Lennard-Jones potential for the fuller- 
ene-nanotube interaction are obtained as the average 
values of the parameters [26] for fullerene-fullerene and 
fullerene-graphene interactions, in accordance with the 
procedure described in [26]. Here, the Lennard-Jones 
potential is used for calculating the potential surface of 
the interaction energy £ w between the fullerene and the 
infinite nanotube, and we believe that this gives ade- 
quate qualitative characteristics of the potential surface 
shape. The cut-off distance, r = r c of the Lennard-Jones 
potential is taken equal to r c = 15 A. For this cut-off 
distance the errors of calculation of the interaction 
energy £ w between the fullerenes and the (8,8) nano- 
tube and the barriers for relative motion and rotation of 
the fullerenes inside the nanotube are less than 0.1%. 
Both the fullerenes and the nanotube are considered to 
be rigid. An account of structure deformation is not 
essential for the shape of the potential surface both for 
the interwall interaction of carbon nanotubes [25,27] 
and the intershell interaction of carbon nanoparticles 
[13,14]. For example, the account of the structure defor- 
mation of the shells of C 60 @C 2 4o nanoparticle gives rise 
to changes of the barriers for relative rotation of 
the shells which are less than 1% [13,14]. It should also 
be noted that the symmetry of interaction energy as a 
function of coordinates describing relative positions of 
interacting objects is determined unambiguously by 
symmetries of the isolated objects and does not change 
if the symmetries of the objects are broken because of 
their interactions. 

The ground state interaction energies between the C 20 
and Fe@C 2 n fullerenes, and the (8,8) nanotube with 
Kekule structure are found to be -1.596 and -1.598 eV, 
respectively. The angles between the C 2 symmetry axes 
of the C 2 o and Fe@C 20 fullerenes and the nanotube axis 
at the ground states are 49.6° and 53.1°, respectively. 
The metastable states with the C 2 symmetry axes of the 
fullerenes perpendicular to the nanotube axis are also 
found for both C 20 and Fe@C 20 - At the metastable 
states, the interaction energies are greater by 4.58 and 
3.04 meV than the ground state energies, for C 20 and 
Fe@C 20 , respectively. 

The potential surfaces of the interaction energy 
between the C 20 fullerene and the nanotube, E w ((p, z), 



as functions of the relative displacement of the fullerene 
along the axis of the nanotube z and the angle of rela- 
tive rotation of the fullerene about the axis of the nano- 
tube (p are presented in Figure 3 for both the considered 
structures of the (8,8) nanotube. In the general case, dif- 
fusion of a fullerene along the nanotube axis is accom- 
panied by rotation of the fullerene. Our calculations 
show that the barriers for rotation of both fullerenes 
about the axes which are perpendicular to the nanotube 
axis lie between 3 and 23 meV for any orientation of 
the fullerene and for both the considered structures of 
the (8,8) nanotube. These barriers are significantly 
greater than the barriers for rotation of the fullerenes 
about the axis of the nanotube (shown in Figure 3 for 
the C 20 fullerene). Therefore, diffusion of the fullerenes 
along the nanotube axis is accompanied only by rotation 
about the nanotube axis. Thus, the minimal barrier A£ d 
for diffusion of the fullerenes along the nanotube axis is 
the barrier between adjacent minima of the potential 
surface E w ((p, z). Figure 3 shows that the shapes of the 
potential surface £w(<P> z)> corresponding to the Kekule 
structure of the (8,8) nanotube and to the structure of 
metallic phase are essentially different. For a case of the 
structure corresponding to the metallic phase of the 
(8,8) nanotube, all the barriers between adjacent minima 
of the potential surface E w ((p, z) are equivalent. In this 
case, the same barrier A£ d = A£ r should be overcome 
for diffusion of the fullerenes along the nanotube axis 
and for rotation of the fullerenes about this axis (see 
Figure 3b). For a case of the Kekule structure, two dif- 
ferent barriers between the adjacent minima exist: the 
barrier A£ d for diffusion of the fullerenes along the 
nanotube axis, and the barrier A£ r to rotation of the 
fullerenes about this axis (see Figure 3a). For the 
Fe@C 2 o fullerene, the shapes of the potential surfaces 
E w (q>, z) are qualitatively the same as for C 20 for both 
the considered structures of the nanotube. The calcu- 
lated values of the barriers A£ d and A£ r are listed in 
Table 1 for both fullerenes. The dependences of the 
interaction energy between the C 2 o and Fe@C 2 n fuller- 
enes and the nanotube on the relative displacement of 
the fullerene along the axis of the nanotube and the 
angle of its relative rotation about this axis, correspond- 
ing to both the considered structures of the (8,8) nano- 
tube are compared as shown in Figure 4. Figure 4 is a 
vivid illustration of the significant changes of the bar- 
riers A£ d and A£ r . The most dramatic change corre- 
sponds to the rotation of the Fe@C 20 fullerene about 
the nanotube axis. Let us discuss the reason of the sig- 
nificant changes of the barriers A£ d and A£ r at the 
Peierls transition by the example of the barrier A£ r for 
rotation of the Fe@C 20 endofullerene about the nano- 
tube axis. Figure 5 presents the dependences E wi ((p) of 
the interaction energies between the endofullerene and 
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Figure 3 The interaction energy £ w (in meV) as a function of the fullerene C 20 displacement z along the nanotube axis and the angle 
q> of the fullerene rotation about the nanotube axis, (a) The (8,8) carbon nanotube with the Kekule structure; (b) the (8,8) carbon nanotube 
with the structure of metallic phase. The energy is given relative to the energy minima. The equipotential lines are drawn at an interval 0.2 meV. 



individual atoms of the nanotube on the angle (p of rota- 
tion of the fullerene. Figure 5 shows that maxima of 
dependences E wi ((p) for individual atoms of the nano- 
tube occur at different angles <p m ; and so the depen- 
dence E w ((p) of total energy on the angle of rotation is 
essentially smoothed. In other words, the barrier A£ r in 
the dependence E w ((p) of the total interaction energy 
between the endofullerene and the nanotube is less by 
an order of magnitude than the barriers A£ rl - in the 
dependences of the interaction energy between the 
endofullerene and only one of the nanotube atoms. 
Thus, the barrier A£ r is very sensitive to the values of 
the barriers A£ ri and angles <p m ;. Therefore, the barrier 
A£ r changes considerably at the found changes of the 
barriers A£ ri and angles <p m , obtained from the Peierls 

Table 1 Calculated characteristics of the dynamical 
behaviour of the C 20 and Fe@C 20 fullerenes inside the 
(8,8) nanotube of different structure 



Nanotube structure 


Kekule structure 


Structure of metallic phase 


Fullerene 


c 20 


Fe@C 20 


c 20 


Fe<aC 20 


Af d (meV) 


1.68 


1.73 


0.87 


0.44 


Af r (meV) 


0.33 


0.05 


0.87 


0.44 


v d (GHz) 


72.2 


54.8 


75.7 


51.3 


v z (THz) 


0.439 


0.144 


0.558 


0.329 


v, (THz) 


1.30 


1.25 


1.36 


1.13 


v y (THz) 


2.15 


2.82 


2.08 


2.14 



A£ d is the barrier for diffusion of the fullerenes along the nanotube axis, A£ r is 
the barrier for rotation of the fullerenes about the nanotube axis, v d is the 
frequency of small relative vibrations of the fullerenes along the nanotube 
axis, v z is the frequency of small relative rotational vibrations of the fullerenes 
about the nanotube axis, v x and vy are the frequencies of small relative 
rotational vibrations of the fullerenes about the two mutually perpendicular 
lateral axes. 



distortions of the nanotube structure and the change of 
the nanotube symmetry at the transition. It should also 
be noted that small barriers to relative motion of 
nanoobjects resulting from the compensation of contri- 
butions of individual atoms to the barriers is a phenom- 
enon well studied by the examples of such systems as 
double-shell carbon nanoparticles [13,14], double-walled 
carbon nanotubes [25,27-32] and a graphene flake in a 
graphite surface [33]. The considerable changes of bar- 
riers for relative rotation of shells at small changes of 
shell structure were found also for double-shell carbon 
nanoparticles [13,14]. 

The frequencies of small vibrations of the fullerenes 
along the nanotube axis (v d ), rotational vibrations about 
the nanotube axis (v z ) and rotational vibrations about 
two mutually perpendicular lateral axes {v x , v y ) are also 
calculated and listed in Table 1. The most remarkable 
change of frequency as a result of the structural phase 
transition corresponds to rotational vibrations of the 
Fe@C2o fullerene about the nanotube axis (this agrees 
with the changes of the barriers). 

Dynamical behaviour of molecules inside 
nanotube 

Let us consider the possible changes of the dynamical 
behaviour of the C 2 o and Fe@C 2u fullerenes inside the 
(8,8) nanotube caused by the structural phase transition. 
The Peierls instability transition temperature T P was 
estimated for the (5,5) nanotube to correspond to tem- 
perature range T P — 1-15 K [4,5,8]. Both barriers A£ d 
and A£ r and the thermal energy k%T P are of the same 
order of magnitude at the structural Peierls phase tran- 
sition (see Table 1). Therefore, dramatic changes of the 
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Figure 4 Interaction energy between the C 20 fullerene or Fe@C 20 endofullerene and (8,8) nanotube Dependences of the interaction 
energy F w between the C 20 fullerene (a, b) or Fe@C 20 endofullerene (c, d) and the (8,8) nanotube on the fullerene displacement z along the 
nanotube axis (a, c) and on the angle tp of the fullerene rotation about the nanotube axis (b, d). Solid lines denote nanotube with the Kekule 
structure, dashed lines denote nanotube with the structure of metallic phase. The energy minimum is positioned at f w = 0, z = 0 and q> = 0. 
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Figure 5 Interaction energy between the Fe@C 20 endofullerene and (8,8) nanotube. Dependences of the interaction energy f w , between 
the Fe@C 20 endofullerene and individual atoms of the (8,8) nanotube on the angle <p of rotation of the endofullerene about the nanotube axis 
are denoted by the thin lines. Dependence of the total interaction energy F w between the Fe@C 20 endofullerene and the (8,8) nanotube on the 
angle <p is denoted by the thick line, (a) The (8,8) carbon nanotube with the Kekule structure; (b) the (8,8) carbon nanotube with the structure of 
metallic phase. All energies are given relative to the energy minima. Only dependences f w , with high values of the barriers A£„ are shown. 
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diffusion and drift over these barriers can take place at 
the Peierls transition for the considered pairs of the 
encapsulated molecules and the nanotube. 

Recently, the diffusion and drift in the periodic poten- 
tial surface of the interaction energy dependent on the 
displacement z along the nanotube axis and the angle (p 
of the rotation about the nanotube axis were considered 
for the thread-like relative motion of walls of double- 
walled carbon nanotubes [28,29]. The expressions for 
the diffusion coefficient and mobility of a movable wall 
were obtained [28,29]. In this study, we consider specific 
cases of diffusion of the molecules along the nanotube 
axis and rotational diffusion about this axis. In this case, 
the expressions for diffusion coefficients, D d and D r , 
mentioned above and corresponding to the diffusion 
along the nanotube axis and rotational diffusion about 
this axis, respectively, take the form: 



■ exp 



hTj' 



exp 



A£ r 



(2) 



where Cl d and O r are the pre-exponential multipliers in 
the Arrhenius formula for the frequency of jumps of the 
molecule between two neighbouring global minima of 
the potential surface E w ((p, z), S d is the distance between 
neighbouring global minima for the motion of the mole- 
cule along the nanotube axis, <5 r is the angle between 
neighbouring global minima corresponding to the mole- 
cule rotation about the nanotube axis and k B is the 
Boltzmann constant. The mobility B d for the motion 
along the axis can be easily obtained from the diffusion 
coefficient D d using the Einstein ratio D d IB d = k^T. 
Figure 3 shows that d d = 0.123 nm and S d = 0.37 nm for 
the (8,8) nanotube with the structure of the metallic 
phase and the Kekule structure, respectively, and S r = 
22.5° for the both structures of this nanotube. 

The value of the pre-exponential multiplier D. in the 
Arrhenius formula is usually considered to be related 
with the frequency v of corresponding vibrations. We 
suppose that the ratio Ci/v remains the same for relative 
motion of different carbon nanoobjects with graphene- 
like structure (nanotube walls and fullerenes). For reor- 
ientation of the fullerenes of the C6o@C24o nanoparticle, 
the frequency multiplier D. was estimated by molecular 
dynamics simulations having the value of 650 ± 350 
GHz [13,14]. We expand the potential surface of the 
intershell interaction energy near the minimum using 
the same empirical potential as in [13,14], and calculate 
the frequencies of small relative librations of the shells. 
The calculated libration frequency has the value v » 50 
GHz, an order of magnitude less than that of the fre- 
quency multiplier O. In the estimations of this study, 
we use the values Cl d ~ 10v d and Q. r ~ 10v z for the pre- 
exponential multipliers. 



The temperature dependencies of the diffusion coeffi- 
cients, D d and D r , estimated using expressions (2) are 
shown in Figure 6. The dependence of the interaction 
energy E w (z) on the displacement of the fullerene C 20 
has two different barriers between the neighbour 
minima for the case of the nanotube with the Kekule 
structure (see Figure 4a). Since the frequency of the 
jumps of the fullerene between the neighbour minima 
exponentially depends on the barrier, the contribution 
of jumps over the lower barrier into the total diffusion 
coefficient is disregarded in these estimations. The tem- 
perature range corresponding to the Peierls transition 
temperature estimates T P - 1-15 K [4,5,8] and the cases 
AEJk^T <1 and AE r /k%T <1, where the Arrhenius for- 
mula is adequate, are considered. (For rotation of the 
Fe@C 2 n inside the nanotube with the Kekule structure 
the Arrhenius formula is not applicable (AE r /k%T >1) at 
this temperature range; this case is considered below.) 
Figure 6 shows that the changes of the diffusion coeffi- 
cients, D d and D r , at the Peierls transition can be of 
orders of magnitude. It is of interest that the diffusion 
coefficient D t for rotational diffusion of the C 20 fullerene 
decreases at the Peierls transition with the increase of 
temperature. 

If a molecule is encapsulated inside a nanotube with- 
out a structural phase transition, the jump rotational 
diffusion takes place at low temperatures, AEJk-gT >1, 
and the free rotation of the molecule occurs at high 
temperature, AE r /k^T <1. Orientational melting (a loss 
of the orientational order with an increase of tempera- 
ture) has a crossover behaviour if the structural phase 
transition is absent. Firstly, orientational melting was 
considered for two-dimensional clusters with shell struc- 
ture [34-37] and later for double-shell carbon nanoparti- 
cles [13,14], double-walled carbon nanotubes [25,38,39] 
and carbon nanotube bundles [40]. In the case where a 
molecule is encapsulated inside a nanotube with a struc- 
tural phase transition and the barrier AE r for rotation of 
the molecule is greater for the high-temperature phase 
than for the low-temperature phase, an inverse orienta- 
tional melting (a loss of the orientational order with a 
decrease of temperature) is possible. In other words, the 
inverse orientational melting takes place if the Peierls 
transition temperature lies in the range AE r i < k%T P 
<A£ r h, where A£ r h and A£ r i are the barriers for mole- 
cule rotation corresponding to high-temperature and 
low-temperature phases of the nanotube, respectively. 
For the considered molecules inside the (8,8) nanotube, 
these temperature ranges are estimated to be 3.8 < T P 
(K) <10 and 0.58 < T P (K) <5.1 for the C 20 and Fe@C 20 
fullerenes, respectively (see Table 1). As these tempera- 
ture ranges are in agreement with the Peierls transition 
temperature estimates T P — 1-15 K [4,5,8], we predict 
that the inverse orientational melting is possible for the 
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Figure 6 Dependences of the diffusion coefficients on the 
reciprocal of temperature 1/7 (a, c) The diffusion coefficient D d ; 

(b) the diffusion coefficient D,. (a, b) The C 20 fullerene; (c) the 
Fe@C 20 endofullerene. Solid lines: the (8,8) carbon nanotube with 
the Kekule structure; dashed lines: the (8,8) carbon nanotube with 
the structure of metallic phase. 



systems considered. The inverse orientational melting 
should be more prominent for the case of the Fe@C 2 o 
fullerene with the greater ratio of the barriers A£ rh /A£ r i. 

Let us discuss the possibility of observing the changes 
of the dynamical behaviour of molecules inside armchair 



carbon nanotubes at the Peierls transition. We believe 
that the most promising method is high-resolution 
transmission electron microscopy. This method was 
used for visualizing dynamics of processes inside nano- 
tubes, such as reactions of fullerene dimerization with 
monitoring of time-dependent changes in the atomic 
positions [41] and rotation of fullerene chains [42]. The 
rotational dynamics of C 6 n fullerenes inside carbon 
nanotube was studied also by analysing the intermediate 
frequency mode lattice vibrations using near-infrared 
Raman spectroscopy [43]. The orientational melting in a 
single nanoparticle may be revealed also by IR or 
Raman study of the temperature dependence of width of 
spectral lines. A specific heat anomaly in multiwalled 
carbon nanotubes may be caused by the orientational 
order-disorder transition [44]. In the case of encapsu- 
lated magnetic molecules (for example, the Fe@C 20 
endofullerene considered above), the study of the tem- 
perature dependence of the electron spin resonance 
spectra could yield information on the molecule rota- 
tional dynamics of these molecules [45]. 

Conclusive remarks 

In this letter, we consider the changes of dynamical 
behaviour of fullerenes encapsulated in armchair carbon 
nanotubes caused by the Peierls transition in the nano- 
tube structure by the example of the C 2 o and Fe@C 20 
fullerenes inside the (8,8) nanotube. We apply the DFT 
approach to calculate the structure of the C 2 n and 
Fe@C 20 fullerenes. The ground state of the (8,8) nano- 
tube is found to be the Kekule structure using the 
method of molecular orbitals. The Lennard-Jones poten- 
tial is used for calculating the barriers for motion of the 
fullerenes along the axis and rotation about the axis of 
the (8,8) nanotube with the Kekule structure and the 
structure with all equal bonds corresponding to low- 
temperature and high-temperature phases, respectively. 
We show that the changes in the coefficients of diffu- 
sion of the fullerenes along the nanotube axis and their 
rotational diffusion at the Peierls transition can be as 
much as several orders of magnitude. The possibility of 
the inverse orientational melting at the Peierls transition 
is predicted. The analogous changes of dynamical beha- 
viour are also possible for other large molecules inside 
armchair nanotubes. We believe that the predicted 
dynamical phenomena can be observed using high- 
resolution transmission electron microscopy, near-infra- 
red Raman spectroscopy, specific heat measurements, 
and by study of electron spin resonance spectra for 
magnetic molecules. 



Abbreviations 

DFT density functional theory. 
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